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ABSTRACT 

In interacting binaries, comparison of a donor star's radial response to mass loss with the response of 
its Roche radius determines whether mass loss persists and, if so, determines the timescale and stability 
of the ensuing evolutionary phase. For giants with deep convective envelopes, the canonical description 
holds that once mass transfer begins it typically proceeds catastrophically on the dynamical timescale, 
as the star cannot lose sufficient heat in order to avoid expansion. However, we demonstrate that the 
local thermal timescale of the envelope's superadiabatic outer surface layer remains comparable to 
that of mass loss in most cases of "dynamical" mass loss. We argue therefore that if mass loss proceeds 
on a timescale longer than this, then even a deep convective envelope will not dramatically expand, 
as the surface layer will have time to relax thermally and reconstitute itself. We demonstrate that in 
general the polytropic approximation gives much too strict a criterion for stability, and discuss the 
dependence of the donor's response on its radius in addition to its core mass. In general, we find 
that the effective response of the donor on rapid timescales cannot be determined accurately without 
detailed evolutionary calculations. 

Subject headings: binaries: close — stars: evolution — stars: mass-loss 



1. INTRODUCTION 

For evolved giant branch stars in binary systems the 
canonical view holds that, should the giant fill its Roche 
lobe (RL), the ensuing mass loss (ML) would in most 
circumstances be dynamically unstable. This is due to 
the donor's deep convective envelope, as it is generally 
held that a convective envelope invar iably expands upon 
ML (|Paczvihski & Sienkiewiczl 1 1972ft . At the same time 
the orbit could either shrink, or expand slower than is 
required to keep the expanding donor inside its RL. Run- 
away ML follows and would then lead to an event typi- 
cally referred to as com mon envelope (C E) evolution, as 
originally conceived by iPaczvnsk i (1976) to explain the 
formation of cataclysmic variables. 

The first quantitative limits o n the stability of ML 
in suc h systems were found by iHiellming fc Webbinkl 
(|1987| ). by modelling the donor star as a condensed poly- 
trope - specifically, an n = 3/2 polytrope with a point 
mass at the center to model the core, implying an isen- 
tropic profile throughout the envelope. The rapid expan- 
sion of these models in response to ML further confirmed 
the susceptibility of giant donors to dynamically unstable 
runaway events. 

The calculations were performed in the so called 
adiabatic regime, which in subsequent studies became 
a standard approach to analyze dynamical instability 
with respect to mass transfer (MT) in binaries. It 
is important here to clarify what is assumed under 
the "adiabatic" approximation: a) a (stellar or poly- 
tropic) model is always in hydrostatic equilibrium (HE); 
b) a model keeps its initial entropy profile as a func- 
tion of mass (jHiellming fc WebbinMll987t iGe et al1[2008l 
l201Qafbh . This means that a characteristic timescale r a( j 
is assumed, which is always shorter than the thermal 



timescale of the star r t h, but longer than the dynamical 
timescale Td yn - Therefore, the obtained model is valid 
only so long as the implied tml satisfies this condition 
in all layers. Models with reduced mass obtained as the 
result of adiabatic ML calculations ought to be driven 
out of thermal equilibrium. 

Consequently, the accuracy of the adiabatic model 
rests on the assumption that the outermost layer of the 
giant, in which convection becomes highly inefficient, is 
negligible in determining MT stability. This can be in- 
terpreted as assuming that such a layer is very quickly 
stripped at the onset of dramatic ML, with the entropy 
profile through all remaining layers unable to adjust in 
order to restore thermal equilibrium. 

In this Letter we argue that in this case the shorter 
thermal timescale of the outer envelope would allow it 
to relax at a rate comparable to or even faster than 
both the expected ML or Tdyn, implying that an adia- 
batic approach (either by considering a polytropic, or 
even a realistic stellar structure) would model the reac- 
tion of the envelope rather poorly. We examine the effect 
of accounting for the thermal relaxation of the outer en- 
velope, and find that in many cases the envelope will 
in fact continue to contract with ML despite the deep 
convective envelope. 

2. RESPONSE OF THE DONOR TO ML 

The formalism for the response to ML of a gi- 
ant branch donor (mode lled as a condensed poly- 
trope) is well-establi shed (jHiellming k Webbinkl 119871 : 
iSoberman et al.lll997l ) . Here we briefly review the canon- 
ical treatment for stability relations related to the stellar 
mass-radius response. In such analysis, stability is de- 
termined from considering small pcrturbative variations 
about the equilibrium state, in which the donor just fills 
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its RL, giving: 



AC 



M SAR 
~R. SM 



(1) 



where C is the response exponent from the mass-radius 
relation (R oc -M^), M and R are the donor mass and 
radius, and Ai? = R — i?RL is the difference between the 
RL radius -Rrl and the donor radius. Depending then 
on the strength of the donor radius' response, ML may 
proceed on one of three timescales: dynamical, thermal, 
or nuclear (t„ uc ). 

For adiabatic ML, a giant star with an isentropic enve- 
lope is modelled as a conde nsed polytrope; using this 
iHiellming fc Webbinkl (|1987l ) found that the adiabatic 
mass-radius exponent of a giant branch star follows to 
good approximation: 
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where m c is the core-mass fraction of the donor. Com- 
parison of Cad with the radius-mass exponent of the RL 
(Crl) then provides the condition on the binary mass ra- 
tio, above which it is believed dynamically unstable MT 
occurs. 

As was mentioned in §1, studies of adiabatic ML thus 
far assume a state of HE in constructing their stel- 
lar models, implicitly taking adiabatic ML to be at a 
timescalc shorter then the thermal timescale of the donor 
but longer then the characteristic dynamical timescale. 
In modelling this process outer mass shells are removed 
while maintaining a constant entropy profile with re- 
spect to the mass coordinate. However if the ML is 
modeled to be adiabatic as described above, such that 
tml = Tad > Tdyn, then the modelled stellar reaction, by 
definition, is not dynamical, but thermal. It is inher- 
ently misleading then, in this type of analysis, to refer to 
the process as dynamical ML, as is typically done. 

In real stars, mass removal from the outer layer of the 
donor will perturb the HE of the star, to which it re- 
sponds on its dynamical timescale. Indeed, as long as 
the ML proceeds on a timescale tml > "?dyn, the star 
and its envelope remain in HE. In this case, the reaction 
of the radius of the star is determined by the reaction on 
the next, longer timescale, presumably thermal, when 
the internal heat is redistributed in the envelope. Deep 
within the envelope, convection is adiabatic; this is very 
efficient and maintains the flat entropy profile (also, its 
value will be constant if tml < ^th)- However, in the 
upper layers of the envelope the density drops off dra- 
matically; convection here is much less efficient and ra- 
diative cooling wins out; the outer reaches of the giant 
become supcradiabatic (V > V a( j). The timescale to re- 
distribute heat in this upper layer could be significantly 
different from that estimated for deep regions based on 
adiabatic convection. It is the relaxation of the outer 
layers which drives MT, giving the relevant timescale at 
high M (jPodsiadlowski et al.ll2002h . 

It is important to emphasize that if the ML occurs ei- 
ther on the dynamical timescale, or on a timescale longer 




Fig. 1. — Entropy (lower panel) and local thermal timescale (up- 
per) in the outer 10% mass of a 1 Mq (R = 30 Rq), 5 Mq 
(R = 50 Rq) and 20 M (R = 950 Rq) giant (black dashed, blue 
dotted and red solid lines respectively). Here, thermal timescale 

was calculated using stellar models as r t ^ i oc = J u(m)dm/ L, 
where u is the specific internal energy. Asterisks indicate the loca- 
tions to the right of (above) which T t ^ i oc < r^ yn . 

than the time necessary for the outer layers to obtain 
their thermal equilibrium, the adiabatic approach is in- 
trinsically invalid. 

Let us now consider realistic stellar models. Within 
real giant branch stars, the entropy profile is marked by a 
steep drop-off in the outermost supcradiabatic "surface" 
of the envelope (see Figure [T]). In prior work concerning 
adiabatic ML it is generally assumed that this layer is 
quickly removed, exposing the convective (dS/dM = 0) 
layer beneath. We note that it is only in this case that the 
donor's structure is well-approximated by a condensed 
polytropic model, in which case on e expects dramatic ex- 
pansion with ML (|Ge et al.ll2010al ). It can be seen as well 
that since in massive giants the superadiabaticity plays a 
role for deep layers as well, the adiabatic approach that 
assumes a flat entropy profile is also invalid. Also, for 
all giants, the local thermal timescale of the supcradia- 
batic layer r t h.ioc is shorter than the donor's dynamical 
timescale: during the thermal relaxation of this layer, 
assumption of HE is not valid a priori. 

For real stars, it is not immediately clear that we can 
safely ignore the response of the surface of the envelope 
to ML. In principle, the rate of ML must surpass the ther- 
mal timescale of the superadiabatic outer envelope in or- 
der for this layer to be effectively removed, and the stan- 
dard picture of adiabatic ML to take hold. Otherwise, 
this layer would have time to thermally relax, reconsti- 
tuting itself faster than it can be torn away. The ther- 
mal timescale of the superadiabatic outer shell, roughly, 
is ~ GMTOghcii/ RL (here m s heii is the mass of the shell 
that has a non-isentropic profile and L is the star's lumi- 
nosity) and the MT timescale as tml = M/M, We then 
have a minimum ML rate that would allow for this outer 
surface of the envelope to be stripped away: 
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where with substantially lower MT rates any response 
of the donor to ML would be poorly modelled using a 
polytropic model. With this limit, a 5 M & giant with 
m shcll - O.OOO1M , L ~ 1OOOL , and R ~ 50i? Q would 
need a minimum ML rate of Af cr ; t ~ 16M Q /yr in order 
to strip the superadiabatic surface layer. For ML rates 
less than this, we cannot assume the envelope's surface 
layer is removed. 

To test the response in real stellar models, we per- 
formed simulations using two different stellar evolution- 
ary codes. The first is a binary stellar-evolu tion code 
STAR S /ev, o rigina l ly developed by Egg leton (|Eggletonl 
fl97Tl 11971 Il97l lEggleton et a TT 1197 3) and then up- 
dated by many other s, as detailed in fPols et alJll995t 
iGlebbeek et al.1 [2008L and references therein.). It is a 
non-Lagrangian code that only calculates stellar model 
in HE, but is capable of treating stars out of thermal 
equilibrium. MT is treated in a parametrized form using 
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where C is 
from 1 to 



a free parameter with values ranging 
10 4 (van der Sluys fc Glebbeek pr i- 
vate communication, see alsolNelson fc Eggleton! (|2001l ); 
lEggleton fe Kiseleva-Eggletonl (|2002[ )). The second is a 
standard Henyey-t ype code, as derived from that origi- 
nally developed in iKippenhahn et all (11967fl: its curren t 
state and input physics are described in llvanoval (|2003 ): 
llvanova fc Taaml ([2004). This code is Lagrangian, and 
is capable of treating stars out of both hydrostatic and 
thermal equilibrium. Here, MT in a binary is calculated 
implicitly keeping the donor always within its RL. 

Using this we can apply a constant rate of ML and ob- 
serve the star's response. As an example, consider a 5Mq 
star on the giant branch with a radius of ~ 50i?©. Apply- 
ing ML rate M = lO _2 M0yr _1 (the canonical adiabatic 
ML regime for this star would be M > 10 -4 M yr _1 ) 
we see that the radius shrinks with ML, despite the deep 
convective envelope of the donor. Note that this result is 
independent of the code itself, as we find the same effect 
using both codes mentioned above. The reason for this 
becomes clear in examining Figure [2] - this is because 
a superadiabatic outer layer remains present even after 
the giant has rapidly lost 0.6M Q . The only difference 
we found is that with the second code m s h is smaller, 
providing greater stability (see Equation [3]) 

Rather than being negligible, understanding the re- 
sponse of this surface layer remains an important factor 
in determining the true response of a giant donor's radius 
to ML. In order to better understand the impact of this, 
we must examine how this impacts MT stability in real 
binaries. 

3. DISCUSSION 

A differential change in the donor radius can be ex- 
pressed as: 
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Fig. 2. — Entropy profile of a 5Mq giant undergoing steady ML 
at a rate of 10~ 2 MQyr~ 1 at the onset (black, solid), and after 
losing ~ 0.6Mq (blue, dashed). 
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where the second term describes the variation in the 
donor's radius with ML over short time intervals, im- 
plying constant entropy and composition, and gives the 
adiabatic response of the donor. The first term can be 
broken up into the normal evolution of the star, and the 
response of the radius during the star's adjustment on its 
thermal and dynamical timescales. The former is negli- 
gible on any timescale considered here. Allowing for the 
relaxation of the outer surface of the envelope, we then 
have an effective mass-radius exponent: 
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For Tth.ioc T ad, the thermal relaxation of the donor's 
superadiabatic surface can dramatically impact the effec- 
tive response of the donor. However, calculations show 
that even for extremely rapid ML (10M Q yr _1 ) this outer 
layer is not fully removed. For MT rates near the crit- 
ical rate estimated by eq. [3] the donor radius expands, 
but not nearly to the dramatic extent we would expect 
in the polytropic approximation: the maximum expan- 
sion is ~ 0.15% of the initial radius (see Figure[3]), while 
in the case where the superadiabatic layer is removed, 
it i s expected to ex pand by ~ 30% (see e.g. Figure 6 
in iGe et all l2010af) - the relative radial expansion in 
MT calculations with real stllar models is 200 times less! 
Well below M cr it, although within what is typically con- 
sidered the adiabatic regime, the donor contracts. 

The consequences for MT stability in a binary are sub- 
stantial. For a 5Mq giant donor with R = 50-R© and a 
core mass of m c = O.856M , the condensed polytropic 
approximation for the donor's response predicts unstable 
MT for any mass ratio q = M/M comp > 0.856 (A/ comp is 
the companion mass). However, with the binary stellar- 
evolution code STARS/ev, the stable MT proceeds for 
mass ratios as high as q = 1.47 (with C = 1), down 
to the nearly complete removal of the original envelope 
(< 1%)Q. 

1 We provide in this stability analysis the results of STARS/ev as 
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Fig. 3. — Mass-radius response for varying ML rates applied to 
an initially 5 Mq giant with radius 50 R@. 

The values given above are not hard limits; using 
STARS/cv, it is difficult to decide in principle when 
any given run has a true run-away, and when the prob- 
lem is numerical. In fact, for 10M Q at 300-R©, we are 
unable to evolve the system through RL overflow even 
below gHj,crit without AR > 10%, which we a ssume 
leads to runaway (see iPodsiadlowski et all 120021 ). The 
parametrization given in eq. [4] is crude at best, and in 
general a better understanding of the degree to which a 
donor may overflow its RL (and the resulting MT rate) 
are needed. We note that for the maximum value of C = 
10000 typically used in STARS (Eggleton 2010, private 
communication), we find g cr it ~ 1.25, although high val- 
ues of C are known to cause the model to crash at the 
onset of MT (van dcr Sluys & Glcbbbcck 2010, private 
communication) . 

Continuing as above, we find a much broader range of 
stability against the onset of a dynamical run-away than 
previously thought, as demonstrated in Table Q] In these 
simulations, we assume MT proceeds in a stable manner 
if the program terminates with either all of the original 
envelope being removed, or at least reaching the end of 
the first episode of MT (for 2.OM < M d < 15My inter - 
mediate mass case B evolution, see IShore et al.lll99l . 
We note as well that under the polytropic approximation 
a binary with q = 1 co uld only avoid a dynamical insta- 
bility for m c > 0.4567 (jHiellming fc WebbinldlT987h : this 
condition does not appear in our calculations (see Fig- 
ure H Table 1). 

Defining £ c ff from simulations is non-trivial (see Fig- 
ure U) : after an initially weaker expansion than antic- 
ipated, C e ff quickly rises before levelling off to a slow 
steady rise through the loss of the first 10% of the donor's 
mass. In Table [U Coff is an averaged value of ( for this 
phase, during which the binary evolves on the quasi- 

it gives least stability in binaries compared to the second code. 



thermal timcscalc (having reached the peak of the Tth 
phase of case B MT, with M ~ 10" 3 M Q /yr). We find 
as well that £ e ff decreases as q increases, therefore the 
value that corresponds to q C rit, reported Ccff.min, is the 
minimum over all stable mass ratios. 

From Table Q] it is evident that any attempt to 
parametrize the adiabatic response in terms of the core 
mass alone will prove unsuccessful. The dependence of 
the donor's response on its radius has been noted in adi- 
abatic models as well (Ge et al, in prep), and suggests 
that the correct response of giant branch donors can only 
be properly computed from detailed stellar evolutionary 
models. At the same time, it is difficult to generalize 
from the result of any such model an absolute criterion 
for MT, as it is difficult to determine when a true run- 
away occurs (see above). In principle the values for g cr j t 
in Table 1 represent only minimum values for the given 
value of C used in eq. 2J 
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Fig. 4. — Initial evolution of the effective radial response of 
the donor (solid lines) and M (dashed lines) during ML, for a 
5Mq-|-5Mq binary with a donor radius of 50Rq (black, highest 
peak), 68-R0 (blue, middle) and 85ii0 (red, lower) with a point 
mass companion. 



Greater stability is found for fully conservative ML 
- no wind losses and all material transferred via RL 
overflow is accreted. Any non-conservation only im- 
proves MT stability through the RL resp onse (Crl) 
ijKalogera fc Webbinkl[T996l lHan et al J 127)021 ) . This is of 
most importance in binaries with a compact accretor, as 
MT there would in the initial stage quickly exceed the 
Eddington limit. If non-conservation of MT is taken into 
account, g C rit can be even higher; for a donor of 5 Mq 
and R ~ 50 Rq we find q cr it ~ 1.61, assuming all mass 
lost from the system leaves with the angular momen- 
tum of the accretor in a bipolar outflow. In a similar 
case, polytropic models would predict t hat only systems 
with q < 1.12 are stable (calculated as in lSoberman et aLl 
Il997f ). Accounting for stellar wind also improves stabil- 
ity further (see Table 1 for 20M Q case, although note 
that at 20Mq the situation is further complicated by 
a much greater depth of the envelope being supcradia- 
batic). Wind losses may also be substantial ly enhanced 
by the companion (|Tout fc Eggleton! 1 19881). rendering 
MT far more stable or even avoidable. lEggletonl (|200l 
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TABLE 1 

Effective critical mass ratio for stable MT 



Donor 






H-W 




Effective 






Ma \Mm\ 

cl L ■ J 




Tllr /A'Ia 




Qcrit 
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9 t 




1 


30 


0.312 


0.158 


0.861 


-0.04 


1.19 
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75 


0.380 


0.284 


0.919 


-0.19 


1.11 
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110 


0.420 


0.368 


0.959 


-0.20 


1.04 




5 


50 


0.171 


-0.065 


0.758 


0.33 


1.47 




5 


68 


0.171 


-0.065 


0.758 


-0.02 


1.11 




5 


85 


0.171 


-0.065 


0.758 


-0.18 


1.02 
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10 


300 


0.223 


0.013 


0.787 


-0.58 


0.64 




10 


350 


0.233 


0.013 


0.796 


-0.35 


0.85 




20 


950 


0.308 


0.151 


0.86 


-0.20 


1.05 




20(18.7) 


950 


0.308 


0.165 


0.86 


-0.19 


1.18 
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our effective response from evolutionary calculations, and the re- 
sulting maximum mass r atios for stable MT. Cad as found from the 
approximation given in ISoberman et all i 19971 ) . Effective values 
for q cr i t and £eff m i n found for C = 1. 



argues that observational evidence suggests post-CE bi- 
naries may only be formed from progenitors where q > 4. 
Even for those systems in which run-away MT (and pre- 
sumably a CE phase) appear to remain inevitable, at 
least 5-10% of the donor's mass may be transferred to 
t he companion prio r to such a phase. 

IGe et all (|2010a[) note that their model for adiabatic 
ML must be understood with the caveat that realisti- 
cally mass would only be lost through a small region in 
the vicinity of the inner Lagrangian point; fully 3 di- 
mensional calculations are needed in order to accurately 
describe the response of the donor. Our studies obey sim- 
ilar limitations, but show the difference in the obtained 
stability of the MT within a similar framework. 

4. CONCLUSIONS 

While adiabatic ML models suggest that it is feasible 
for rapid ML to drive dramatic expansion of the envelope, 



this rests on the assumption that the supcradiabatic sur- 
face can be removed completely. 
We conclude that: 

thermal relaxation of the supcradiabatic surface is 
reflected in the effective mass-radius response of the 
donor and affects the resulting stability strongly; 



no parametrization can be done based only on q 
and/or on the ratio of the core mass to the donor 
mass; 



the response of the donor necessarily evolves with 
ML, rather than remaining static. 

In principle one can only make strong conclusions for 
the response of any given donor on a case by case basis, 
using a detailed stellar evolutionary code. At the same 
time, one must acknowledge that attempting to deter- 
mine Q cr it for any given donor/core mass using such a 
code is inherently model-dependent. 

Certainly ML from giant branch donors is more stable 
than previously thought from simplified adiabatic mod- 
els. It is known that consideration of non-conservative 
MT shows stability is likely greater still, in particular 
for compact accretors in which MEdd is quickly reached. 
However, to make any strong claims regarding the con- 
sequences for the formation rates of any given post-CE 
binary would be rather speculative without further pop- 
ulation studies. 
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